Evolution of Nuclear Shell Structure due to the Pion Exchange Potential 



O 

o 

(N 

(D 

OO 
(N 



> 

o 

O 

\o 

O 



X 



WenHui Long,^' ^'0 Hiroyuki Sagawa,^ Jie Meng,^ and Nguyen Van Giai'^ 

'School of Physics, Peking University, 100871 Beijing, China 
^ Center for Mathematical Sciences, University of Aizu, Aizu-Wakamatsu, 965-8580 Fukushima, Japan 
"^Institut de Physique Nucleaire, CNRS-IN2P3, Universite Paris-Sud, 91406 Orsay, France 



The evolution of nuclear shell structure is investigated for the first time within density-dependent 
relativistic Hartree-Fock theory and the role of 7r-exchange potential is studied in detail. The energy 
differences between the neutron orbits {i'l/!.9/2, ^^1*13/2} in the N = 82 isotones and between the 
proton ones {'iTlg7/2, 7i"l/iii/2} in the Z = 50 isotopes are extracted as a function of neutron excess 
N — Z. A kink around Z = 58 for the N = 82 isotones is found as an effect resulting from pion 
correlations. It is shown that the inclusion of 7r-coupling plays a central role to provide realistic 
isospin dependence of the energy differences. In particular, the tensor part of the 7r-coupling has an 
important effect on the characteristic isospin dependence observed in recent experiments. 
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The nucleon-nucleon bare interaction at low and 
medium energies is originally induced by the meson 
exchange processes as predicted by Yukawa 0. In 
the nuclear medium, it is strongly renormalized by 
medium effects which lead to the effective nucleon- 
nucleon interaction. A large part of the present 
understanding of nuclear structure is based on self- 
consistent mean field descriptions making use of ef- 
fective interactions directly parametrized so as to re- 
produce selected nuclear properties. The two most 
successful categories are the non-relativistic Hartree- 
Fock approaches |^ ^ [3 and the Relativistic Mean 
Field (RMF) approaches [i, 0, Q ■ Using these mod- 
els, many nuclear structure properties are calculated 
in the whole region of the nuclear chart with these ef- 
fective interactions. During the past decade, great 
successes have been achieved by the RMF theory 
not only in stable nuclei but also in exotic regions 
|ilISIllllll3,,I4lIllIllI3- Of special interest 
is the fact that the RMF model provides a natural 
mechanism for explaining the spin-orbit splittings of 
single-particle levels. This feature becomes even more 
of central importance with the experimental observa- 
tion that nuclei near drip lines undergo modifications 
of their shell structure, where the spin-orbit potential 
must play an essential role. 

One of the basic open problems is the role of 
one-pion exchange process, which is known to play 
a fundamental role in the meson-exchange interac- 
tion. However, the RMF model is not the appropri- 
ate framework to study pion-related processes because 
it is essentially a Hartree approximation where the 
Fock (exchange) contributions are altogether dropped, 
while the Hartree (direct) contributions of pions are 
zero due to the parity conservation in spherical and 
axially deformed nuclei. Recent progress in the rela- 
tivistic Hartree-Fock description of nuclear structure, 
namely the density-dependent relativistic Hartree- 
Fock (DDRHF) approach [H has brought a new in- 
sight to consider this problem. Within the DDRHF 



theory, the effective meson-nucleon coupling strengths 
including the one-pion exchange are determined in a 
similar way to the RMF model and a quantitative de- 
scription of the nuclear structure properties can be 
successfully achieved ^jll^. Thus, DDRHF opens 
the possibility to investigate the role of one-pion ex- 
change processes in nuclear structure problems within 
the framework of a relativistic mean field theory. 

In this work, we study the role of pion-exchange 
processes on single-particle spectra by concentrating 
on specific cases. In the recent paper of Schiffer et 
al. [21| , it was shown that a set of states outside the 
proton core Z = 50 and the neutron core = 82 
may provide a unique information to determine the 
evolution of the nuclear shell structure. This is why 
we choose here to discuss the pion effects taking these 
nuclei as an illustrative example. 

The one-pion exchange process with pseudo-vector 
coupling contains two types of contributions, the cen- 
tral coupling and the non-central tensor coupling |22| . 
Recently, the tensor type force was shown to have a 
distinct effect on the evolution of the nuclear shell 
structure in a non-relativistic Hartree-Fock model 
[23I I . In this work, we will study the behavior of single- 
particle energies of the states {i/l/ig/2, 1^^113/2} {1^ de- 
notes neutron states) in the iV = 82 isotones and the 
states {7rl(77/g, Trl/in/a} (tt denotes the proton states) 
in the Z = 50 isotopes within the DDRHF theory. 
The effect of the 7r-coupling, especially the contribu- 
tion of its tensor component on the isospin dependence 
of the shell evolution will be analyzed in detail. 

In a non-relativistic reduction, the one-pion ex- 
change potential can be divided into two parts, V^{q) 
and Vj{q): 
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The Fourier-transform of the central part gives a re- 
pulsive contact interaction and an attractive Yukawa 
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In general, the repulsive contact part is strongly hin- 
dered by the Pauli principle between two nucleons. 
Moreover, in DDRHF the u;-nieson exchange will 
strongly suppress the repulsive contact part so that 
it can be discarded on a solid basis The central 
part of the 7r-induccd potential is thus given by 

^\Y1 Ca4'^'/^'5 / dridr2fa{ri)fi3{r2) 

afi-rS (3) 

X K«^](l,2)/7('^2)/5(ri), 

where Ci and cj correspond to nucleon annihilation 
and creation operators, and the interaction vertex in- 
cluding the propagator is 



(1,2) 



g-m,|ri-r2| 
127r |ri - r2| 



(4) 

The tensor 7r-coupling Vj can be obtained similarly. 

In the rest of this paper, the DDRHF calculations 
are performed with the parameter sets PKOl, PK02, 
PK03 given in Ref. , and they are compared with 
RMF (Hartrec) results calculated with the parameter 
set PKDD O. It should be noticed that PK03 has 
a stronger 7r-coupling than PKOl, while PK02 and 
PKDD have no 7r-coupling in the Lagrangian. The 
pairing correlations are treated by the BCS method 
with a density-dependent, zero-range pairing force 
[2^ 1 . The nuclei studied are located along the isotones 
with iV = 82 from ^^^Sn to ^^^W and the Sn isotopes 
with Z = 50 from ^°'^Sn to ^^^Sn. 

The energy differences AE{Z) between two unoc- 
cupied neutron states |j/l/ig/2, ^^1*13/2} of the = 82 
isotones are shown in the first panel of Fig. [Jwhilc the 
values AE{N) between two unoccupied proton states 
{7rl(77/2, 7rl/iii/2} of the Z ~ 50 isotopes arc shown 
in the second panel. The results of pion-dependent 
models PKOl and PK03 show kinks at Z = 58 
{N - Z = 24, i.e.,"OCe) for the iV = 82 isotones. 
In contrast, the pion-independent models PK02 and 
PKDD predict a monotonous decrease with respect to 
N — Z. The predictions of PKOl and PK03 give more 
realistic isospin dependence, qualitatively and quanti- 
tatively, than those of PK02 and PKDD as compared 
with the experimental data, although substantial dis- 
crepancies still exist. A similar isospin dependence 
is also observed for the proton states of the Sn iso- 
topes in the second panel of Fig. [1] the inclusion 
of TT-coupling in RHF (PKOl and PK03) brings also 
a significant improvement compared with the results 
calculated without 7r-coupling (PK02 and PKDD). 
Around ^^^Sn there exists a small kink in the data 
and only PK03 gives such a trend, which might be 
due to its stronger 7r-coupHng than PKOl. 

The single-particle energies contain contributions 
from the different mesons. The pion contributions can 
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FIG. 1: The energy differences AE{Z) = E^u^^^^ - 
£■^1/19/2 in A'' = 82 isotones (first panel), and AE{N) = 
-^t1''ii/2 ~ ^Ti9~/, in Z = 50 isotopes (second panel) as a 
function of neutron excess. The DDRHF results are due 
to PKOl, PK02, PK03, while the RMF results are ob- 
tained with PKDD. The experimental data are extracted 
from Ref. El. 



in turn be separated into a central and a tensor com- 
ponent. In Fig. 12 arc shown the 7r-coupling contribu- 
tions to the AE{Z) in the = 82 isotones calculated 
with PKOl. Comparing to the first panel of Fig. U 
one can see that the characteristic isospin dependence 
of AE, namely the kink a,t Z = 58, is mainly due to 
the pion coupling and particularly to its tensor part. 
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FIG. 2: Central and tensor pion contributions to the en- 
ergy differences AE{Z) = £'1^1113/2 ~ Evihg/2 'iue to the 
TT-coupfing, as a function oi N — Z for the N = 82 isotones. 
The results are calculated by DDRHF with PKOl. 
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FIG. 3: The contribution of 7r-coupling to the energy dif- 
ference AE = -Biyiii3^2 — Eyihg/2 ^ ^ function of N — Z 
along the isotonic chain A*' — 82. In the first panel, it is 
separated into three parts: Neutron Core (A'^ = 82), Pro- 
ton Core {Z = 50) and Valence Protons. The contribu- 
tions of the tensor 7r-coupling from 4 valence proton states 
are shown in the second panel. The results correspond to 
DDRHF with PKOl. 



We can deepen further our understanding by ex- 
amining the configuration dependence of the tt- 
coupling contributions. From ^^^Sn to ^^^W the 
added protons will gradually occupy the valence or- 
bits 7Tlgy/2,n2d^/2,n2d3/2 and -KXh^xji- In Fig. Owe 
separate the contributions of 7r-coupling to AE{Z) 
into three components from the neutron core orbits 
(iV = 82), the proton core orbits {Z = 50) and the va- 
lence proton orbits, respectively. From the first panel 
of Fig. [31 it is seen that the valence protons play 
the most important role to give the isospin depen- 
dence shown in Fig. [J and Fig. |5J In the second 
panel of Fig. the contributions of 4 valence proton 
orbits to AE{Z) arc separately plotted to examine 
the configuration dependence of the tensor 7r-coupling. 
As the proton number increases from Z = 50 to 58 
(TV — Z = 32 to 24), the 7TlgT/2 orbit is gradually oc- 
cupied. Then, from Z = 60 to 64 (iV - Z = 22 to 18), 
the protons fill the 7r2(i5/2 orbit and finally the 7rl/iii/2 
orbit is occupied in the nuclei with Z = 66 ~ 74 
{N — Z = 16 to 8). It is found from the second panel 
of Fig. Othat the magnitude of the tensor correlation 
is very sensitive to the occupation probabilities of the 
valence orbits, especially, those of high-j orbits. 

To have a better understanding about the contribu- 
tions of these valence proton states, we calculate the 



interaction matrix elements of 7r-coupling, separating 

the central and tensor parts, vj^' '^^ and V^^' Fig. 
13] shows the central (first panel) and tensor (second 
panel) matrix elements between the neutron states, 
1^1/19/2 or i^lii3/2 and the different valence proton 
states. The central 7r-coupling is attractive and it 
shows a simple monotonous isospin dependence for 
both iylhg^2 a-nd i'lii^/2- In contrast, the tensor part 
shows different isospin behaviours depending on the 
configurations. Namely, the two-body matrix ele- 
ments between vlhg/2 and 7Tlgif2 or 712(1^/2 arc re- 
pulsive whereas those between vlhg^2 and 7rl/i]^]^/2 or 
7r2d5/2 arc attractive. For i^lii3/2, they have com- 
pletely opposite signs as seen in the second panel of 
Fig.EI 
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FIG. 4: The central (first panel) and tensor (second panel) 
parts of the interaction matrix elements between neutron 
states (i^l/ig/2 (filled symbols) and^lji3/2 (open symbols)) 
and valence proton states {iTlgT/2 (squares), 7r2d5/2 (cir- 
cles), 7r2ci3/2 (up-triangles) and nlhii/2 (down-triangles)) 
induced by the 7r-coupling as function of A'^ — Z in the 
N = 82 isotones. The results correspond to DDRHF with 
PKOl. 

The first panel of Fig. |3] shows that the central tt- 
coupling is always attractive. As a result, it provides 
the simple monotonous behaviour and weak isospin 
dependence (see Fig. [SJ. At first sight, the tensor tt- 
coupling looks more complicated. However, one can 
find the following regularities: 1) the isospin depen- 
dence of the matrix elements is determined by the oc- 
cupation probabilities of the valence proton orbits. 2) 
the tensor 7r-coupling is attractive for (j< , j> ) config- 
urations, while it is repulsive for (j>, j>) and {.]<,.]<) 
configurations. 3) the magnitude of the tensor intcrac- 
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tion matrix element is proportional to the degeneracy 
(2j + 1) of the orbit. These properties of the tensor 
force in DDRHF with PKOl and PK03 can provide a 
more reasonable isospin dependence of the energy dif- 
ference between 1^1/19/2 and z/1ji3/2. This pecuhar fea- 
ture of the tensor interaction was already pointed out 
in the dcutcron binding mechanism in Ref . 1251 . a nd 
also in the non-relativistic HF calculations |23, |2J| . 

In this work, the evolutions of nuclear shell struc- 
ture outside the = 82 neutron and Z = 50 pro- 
ton closed shells are investigated for the first time 
within the density-dependent relativistic Hartree-Fock 
(DDRHF) theory and the role of the one pion- 
exchange potential is studied in detail. The en- 
ergy differences AE{N — Z) between the states 
{1^1/19/2,1^1113/2} in the = 82 isotones, and those 
between {7rl(77/2, 7rl/iii/2} in the Z = 50 isotopes are 
extracted as a function of neutron excess to explore 
the shell evolution of these exotic nuclei. It is found 
that a better isospin dependence of AE{N — Z) is ob- 
tained if one includes a pion-coupling term in the ef- 
fective Lagrangian. This pion coupling gives the char- 
acteristic isospin dependence observed in recent ex- 
periments. Furthermore, the tensor 7r-coupling with 
the valence protons provides the important physical 
mechanism to give a proper description of the isospin 
dependence. Indeed, the strong configuration depen- 
dence of the tensor part of the 7r-coupling plays a cru- 
cial role to give the kink of the AE{N — Z) observed 



in the experimental systematics of the shell evolution. 

Compared with the experimental results, however, 
there still exist some substantial discrepancies. Since 
the studied states are unoccupied orbitals, the dy- 
namical coupling to the core vibrations which is not 
considered here, might be one missing physical mecha- 
nism to cure the existing disagreement. In the present 
DDRHF the 7r-nucleon effective coupling, determined 
by fitting the empirical data of the ground state prop- 
erties of selected nuclei, is density-dependent and its 
strength turns out to be fairly weak in the medium 
1L9| . Another promising meson coupling for the shell 
evolution is the Lorcntz p-tensor coupling which is 
not included yet. As a perspective, the inclusion of 
dynamical coupling to RPA vibrations and the treat- 
ment of p-tensor coupling might play an additional 
important role to cure the existing discrepancy. 
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